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We investigate effects of Coulomb interaction and hopping transport in the insulator phase of 
granular metals and quantum dot arrays. We consider a spatially periodic as well as an irregular 
array, including disorder in a form of a random on-site electrostatic potential. We study the Mott 
transition between the insulating and metallic states in the regular system and find the dependence 
of the Mott gap upon the intergranular coupling. The conductivity of a strictly periodic array has 
an activation form with the Mott gap as an activation energy. Considering irregular systems we 
concentrate on the transport properties in the dielectric, low coupling limit and derive the Efros- 
Shklovskii law for hopping conductivity. In the irregular arrays electrostatic disorder results in the 
finite density of states on the Fermi level giving rise to the variable range hopping mechanism. We 
develop a theory of tunneling through a chain of grains and discuss in detail both elastic and inelastic 
cotunneling mechanisms; the former dominates at very low temperatures and/or very low applied 
electric fields, while the inelastic mechanism controls tunneling at high temperature/fields. Our 
results are obtained within the framework of the new technique based on the mapping of quantum 
electronic problem onto the classical gas of Coulomb charges. The processes of quantum tunnelling 
of real electrons are represented in this technique as trajectories (world lines) of charged classical 
particles in d+ 1 dimensions. The Mott gap is related to the dielectric susceptibility of the Coulomb 
gas in the direction of the imaginary time axis. 

PACS numbers: 73.23Hk, 73.22Lp, 71.30.-hh 



I. INTRODUCTION 

Electron transport in granular metals is a subject of in- 
tense theoretical and experimental research i&SdjtZ&Z&L. 
The interest is motivated by the fact that granular met- 
als represent a unique tunable system that captures all 
the essential features of general disordered systems and 
allows for the study the effects of interference of strong 
Coulomb interaction and disorder. In particular, depend- 
ing on the intergranular conductance, granular metals ex- 
hibit the wide spectrum of transport behaviors ranging 
from the hopping conductivity to those typical to disor- 
dered Fermi liquids. 

At present the metallic regime corresponding to strong 
intergranular coupling where the Coulomb repulsion is 
substantially screened is understood fairly well 3,5 . The 
key energy scale in this regime is the inverse escape rate 
from the grain T = g5, where S is the mean level spacing 
in the individual granule, and g is the intergranular con- 
ductance. Generally, on the energy scales larger than T 
the granularity is essential, while at lower energy scales 
the system resembles an ordinary disordered metal with 
an effective diffusion coefficient determined by the cou- 
pling between the grains. As a result at temperatures 
larger than F the conductivity exhibits the logarithmic 
temperature dependence in all dimensions, while the low 
temperature phase at T < T, 3 can be viewed as the disor- 
dered Fermi liquid 5,6 with conductivity described by the 
Altshuler-Aronov— and localization type corrections 

The dielectric regime where the conductivity is known 
to depend on the temperature asi^: 

a = a exp[-(T o /T) 1 / 2 ], (1) 

with ctq being the high temperature conductivity and T D 



being the characteristic energy scale defined below, has 
been long posing a puzzle. This form of the observed 
behavior suggested the primary role of Coulomb corre- 
lations and several models were proposed, yet the prob- 
lem of conduction in granular metals remained unsolved, 
see Ref. [l3| for a review. Building on the constructive 
critique of Ref. , the important development towards 
understanding the observed behavior JIJ in granular met- 
als was made in Ref. Q, where the role of electrostatic 
disorder in formation of random potentials on individual 
granules was fully recognized and explored. This sug- 
gested that the dependence can indeed be viewed 
as Mott-Efros-Shklovskii variable range hopping in the 
presence of the Coulomb gapii The remaining problem 
however was that in order to realize variable range hop- 
ping conductivity an electron have to tunnel through the 
sequence of granules in one hop and no mechanism that 
could provide such a tunneling process was described. 

In the Letter 15] we have reported the mechanism, 
the multiple coherent co-tunneling, via which the corre- 
lated hopping in a Coulomb gap occurs in granular con- 
ductors. We developed the technique allowing rigorous 
consideration of Coulomb effects putting thus the theory 
of hopping transport in granular conductors on a firm 
quantitative theoretical basis. In the present work we 
investigate in detail the effects of the Coulomb interac- 
tion on transport properties of granular metallic arrays 
concentrating on the regime of low and intermediate in- 
tergranular coupling. We discuss both the model of a 
periodic granular array as well as the irregular model 
that includes the random on-site electrostatic potential. 
We begin with the periodic model and show that the 
regular array of metallic granules is a Mott insulator at 
zero temperature, provided the tunneling conductance is 



2 



less than the critical value g c = (1/7tz)]il(Eq/5) where 
z is the coordination number ( number of neighbors for 
a site on the array ), 5 is the mean energy level spac- 
ing and Eq is the Coulomb charging energy in a single 
grain. The Mott gap, the characteristic property of the 
insulating state, decreases with the increase of the inter- 
granular conductance and disappears when g riches the 
critical value g c . 

We develop an approach to the insulating state of gran- 
ular conductors based on the mapping the original quan- 
tum model describing electrons in the granular system 
on the classical gas of the interacting Coulomb charges. 
More specifically, we first absorb the Coulomb interac- 
tion into the phase field, similarly to Ref. |l6( and next 
map the quantum problem that involves functional in- 
tegration over the phase fields onto the classical model 
that has an extra time dimension - in an analogy with 
the method of Ref. The classical charges of the 

emerging effective model are connected via the electron 
world lines representing graphically tunneling processes 
and form necklace-like loops. In the insulating state only 
short loops are present; this means that electrons and 
holes form bound states and only virtual tunnellings over 
the short distances are allowed. In the metallic phase the 
loops of the infinite size are present meaning that elec- 
trons and holes are not bound together. We show that 
the Mott gap in the insulating state is given by the time- 
axis component of the inverse dielectric susceptibility of 
the classical charged gas. 

Mott transition in a periodic granular array resem- 
bles the Mott transition in the Hubbard model 18 . How- 
ever, there is an important difference between the two. 
Namely, in granular metals, even in the case of the pe- 
riodic granular samples, the electron motion inside the 
grains as well as on the scales exceeding the intergranular 
distance is diffusive, while the Hubbard model assumes 
that electrons move through a periodic lattice that allows 
to label their states by quasimomenta. The model of a 
granular metal also has an extra physical parameter (5, 
the mean energy levels spacing in a single grain, that has 
no analog in the Hubbard model. Yet, the physics of the 
transition is similar in both cases, and one may expect 
that the methods developed for the study of the Mott 
transition in the system of strongly interacting electrons 
on the lattice (such as dynamical mean field theory 18 , for 
example) can be applied to the granular metals as well. 
At present, however, the description of the Mott gap in 
granular arrays in the very vicinity of the Mott transition 
is not available, in particular the order of the transition 
remains an unresolved problem. 

The spatially regular granular arrays, including tech- 
nologically important nanocrystals of quantum dots, are 
now the subject of the extensive experimental study. One 
might expect that periodic arrays would show activa- 
tion transport, characteristic for the systems with the 
hard gap in the electron spectrum. However, even seem- 
ingly perfectly periodic arrays2ii2iSi do not follow the 
expected activation conductivity behavior. Instead, the 



experimentally observed conductivity follows the Efros- 
Shklovskii (ES) law expressed by Eq. QJ- This behavior 
indicates that electrostatic disorder plays a primary role 
in these systems. 

Dependence (QJ cannot be observed in the strictly peri- 
odic structures without disorder—, where the hard gap in 
the excitation spectrum results in the activation behavior 
of the conductivity. The two necessary ingredients for the 
variable-range- hopping type(VRH) conductivity are: (i) 
the finite density of states at energies close to Fermi level 
and (ii) the exponentially falling probability of tunneling 
between the states close to the Fermi level. Furthermore, 
the transformation of the usual Mott's VRH dependence 
into the Efros-Shklovskii law reflects the presence of the 
so-called soft gap in the electron spectrum resulting from 
the strong long range Coulomb correlations. The prob- 
lem of hopping transport in granular conductors is thus 
two-fold: (i) to explain the origin of the finite density of 
states near the Fermi-level and the role of Coulomb cor- 
relations and (ii) to elucidate the mechanism of tunneling 
through the dense array of metallic grains. 

Earlier attempts to explain the behavior JJJ were based 
on the assumption of the large dispersion in the granule's 
size— Moreover, one more strong assumption was nec- 
essary in order to arrive at formula QJ that the distance 
between the grains and the grain sizes are not statisti- 
cally independent but rather are linearly coupled with 
each other, and this hypothesis contradicted to exper- 
imental findings^. Furthermore, recent observations of 
ES law were made in the arrays of quantum dotsiS and ar- 
tificially manufactured granular systems 20 where the size 
of granules/periodicity in the dots arrangement was con- 
trolled within the few percent accuracy. Thus the basic 
assumption about the large size/distance dispersion does 
not apply to all these systems, yet the dependence Q 
is observed. Moreover, as was pointed out in Refs. pll3j 
the capacitance disorder alone can never lift the Coulomb 
blockade in a single grain completely, and thus this mech- 
anism for sure cannot explain the finite density of states 
at the Fermi level, and as a consequence it cannot explain 
dependence at low temperatures. 

To address this problem properly, one notices that the 
insulating matrix in granular conductors that is typically 
formed by the amorphous oxide, may contain a deep tail 
of localized states due to carrier traps. The traps with 
energies lower than the Fermi level are charged and in- 
duce the potential of the order of e 2 / nr on the closest 
granule, where n is the dielectric constant of the insula- 
tor and r is the distance from the granule to the trap. 
This compares with the Coulomb blockade energies due 
to charging metallic granules during the transport pro- 
cess. This mechanism was considered in Ref. 0. If 
granules are very small, surface effects also contribute to 
random changes in associated Coulomb energies. And, 
finally, talking about the 2D granular arrays and/or ar- 
rays of quantum dots, one expects that the source of the 
induced random potential are imperfections of and the 
charged defects in the substrate. Thus in granular con- 
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ductors the role of the finite density of impurity states 
is taken up by random shifts in the granule Fermi levels 
due to all the above reasons. We describe these shifts 
as the external random potential /i;, where i is the grain 
index, applied to each site of a granular system. 

Tunneling over the distances well exceeding the aver- 
age granule size in a dense granular array is a funda- 
mental physical problem. The virtual electron tunneling 
through a single granule (quantum dot) was considered 
in Ref. [2^ where two different mechanisms for a charge 
transport through a single quantum dot in the Coulomb 
blockade regime were identified. The first mechanism - 
the so called elastic cotunneling - transfers the charge via 
tunneling of an electron through an intermediate virtual 
state on the dot such that the electron comes out of the 
dot with the same energy it came in. In the second mech- 
anism - inelastic cotunneling - an electron that comes out 
of the dot has a different energy from the one that came 
in. After such a tunneling process an electron-hole excita- 
tion is left in the dot which absorbs the energy difference. 
Note that both these processes go via the classically inac- 
cessible intermediate states, i.e. both mechanisms occur 
as the charge transfer via a virtual state. Inelastic cotun- 
neling dominates the elastic one in the region where the 
temperature or the applied voltage is larger than ^EqS. 

Coulomb charge mapping technique, developed in our 
work, allows for a simple and transparent analytic de- 
scription of both elastic and inelastic cotunneling pro- 
cesses offering thus a ground for a comprehensive theory 
of the hopping transport in granular metals. The advan- 
tage of our approach is that the higher order tunnelling 
processes can be included in a straightforward way along 
with the nearest neighbor hoppings. The obtained re- 
sults for the hopping conductivity dependencies as well 
as results concerning regular arrays are summarized in 
the section below. 

Our paper is organized as follows: We summarize our 
results in Sec. [H] and introduce our model in details in 
Sec. IIIII In Sec. |^ we find the corrections to the Mott 
gap due to the weak intergranular coupling. In Sec. |3 
we introduce the charge representation for the model of 
a granular metal and discuss the physical quantities in 
this representation in Sec. IVII The Mott gap for the pe- 
riodic sample in the case of moderately large tunneling 
conductance is found in Sec. IVIII In Sec. IVIIII we gen- 
eralize our technique to include the high order tunneling 
processes and apply this technique to find the hopping 
conductivity in Sec. lIXI and to study the Mott transition 
in Sec. |3 Our conclusions are presented in Sec. IXII while 
some technical details are relegated into Appendix. 



II. SUMMARY OF THE RESULTS AND 
DISCUSSION 

To treat the coherent multi-granule tunneling and to 
describe the insulating state we have developed an ap- 
proach based on the mapping of the quantum model of 



granular metal onto the classical gas of Coulomb charges 
in d+l dimensions. The Coulomb charges of the effective 
classical model are connected by electron world lines and 
form necklace-like loops. First we apply our technique to 
calculate the Mott gap Am in the electron spectrum of a 
periodic granular array and show that it is simply related 
to the dielectric constant n r of the classical charge gas in 
the extra time direction. In the insulating regime where 
virtual tunnellings on large distances are suppressed, the 
charge model approach is reduced to consideration of the 
short loops, representing the nearest neighbor hopping. 
Close to the metal-insulator transition the larger loops 
that describe tunneling on larger distances become im- 
portant, and in the metallic regime the loops of an infinite 
size appear. In this picture, the Mott transition corre- 
sponds to the transition between the state where large 
loops are energetically suppressed and the state where 
loops of an infinite size do proliferate. It is important that 
despite the large loops representing tunneling over large 
distances are energetically suppressed and do not con- 
tribute to thermodynamics of the insulating state, they 
are essential for transport and govern hopping conduc- 
tivity in the arrays subject to electrostatic disorder. 

We will discuss the physics of the Mott insulator phase 
and describe the metal to insulator transition in the con- 
text of the periodic granular model. The subsequent 
derivation of hopping conductivity in the weakly coupled 
arrays subject to external electrostatic disorder does not 
require special assumptions about the spatial periodicity 
and holds both for regular and irregular arrangements. 



A. Periodic granular arrays 

An array of weakly coupled metallic grains, g <C 1, is 
the charge insulator at zero temperature. The electron 
transport is blocked due to the presence of the Coulomb 
gap Am in the electron excitation spectrum. In the 
limit of vanishingly small intergranular tunneling conduc- 
tances this gap is given by the Coulomb charging energy 
of a single grain Am = Eft. At finite tunneling conduc- 
tance the Mott gap Am is reduced due to the virtual 
electron tunneling to neighboring grains. Suppression of 
the Mott gap Am due to such virtual processes in the 
limit of weak coupling can be obtained from the pertur- 
bation theory and is given by 

A M = - — 9 Eeh In 2, gz4Zl, (2) 

7T 

where z is the coordination number and E e h is the en- 
ergy to create an electron-hole excitation in the system 
by removing an electron from a given site and placing 
it on one of the neighboring sites. In the case of diag- 
onal Coulomb interaction E^ — Eq Sij the energy of an 
electron-hole excitation E e h is simply twice the Coulomb 
charging energy Eq. From Eq. J2J) it follows that the 
Mott gap Am is reduced significantly at the values of 
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FIG. 1: Schematic dependence of the Mott gap on the tun- 
neling conductance g for a periodic granular array of the di- 
mensionality d > 2. At zero tunneling conductance the Mott 
gap is given by the charging energy in a single grain and it 
suppressed exponentially with g for gz > 1. The Mott tran- 
sition takes place at Am ~ g5, where 5 is the mean energy 
level spacing in a single grain. The possible continuous tran- 
sition that cannot be ruled out at present time is shown by 
the dotted line. 

the tunneling conductance gz ~ 1, where at the same 
time the perturbation theory becomes unapplicable. 

In order to find the Mott gap Am at larger tunneling 
conductances, we make use of mapping of the original 
model onto the model of classical charges and construct 
a mean field theory considering the coordination number 
z as a large parameter, z ^> 1. Within this theory we 
obtain that the Mott gap Am is exponentially suppressed 
with the tunneling conductance 

A M = cgEZe-™9, A M > gS gz>l, (3) 

where c is the number of the order of unity. When de- 
riving equation (JHJ) , the nearest neighbor virtual electron 
tunneling only was taken into account. In the related 
effective classical model this corresponds to accounting 
for the short quadruple loops only. For this reason the 
result © holds only as long as the Mott gap Am is 
larger than the inverse electron escape time from the 
grain T = gS. Solving the inequality Am > gd together 
with the expression for the Mott gap (0) within the log- 
arithmic accuracy we obtain that the later can be used 
as long as the tunneling conductance g is less than the 
critical conductance given by the expression 

g c = — HE C /S). (4) 

Z7T 

In the regime where the Mott gap Am becomes of the or- 
der of the inverse escape time of an electron from a grain 
gS the large loops become important showing tendency 
to the formation of the metallic state. This means that 
the value (|3J) determines the boundary of stability of the 
insulating state. At the same time the critical tunnel- 
ing conductance Q coincides with the boundary of the 
stability of the metallic phase obtained in Ref. [5| mean- 
ing that the critical value (@J) indeed marks the boundary 



between the metallic and insulating phases at zero tem- 
perature. 

Considerations leading to the result Q for the crit- 
ical tunneling conductance g c are based on the investi- 
gation of the behavior of the system from the insulat- 
ing and metallic sides implemented in the present paper 
and Ref. p] respectively. However, both these considera- 
tions do not apply in the immediate vicinity of the phase 
transition. At present the quantitative description of the 
Mott transition is not available, in particular, the order 
of the transition remains an open question. Two possible 
options are illustrated on Fig[I] 

Presented description of the Mott transition may not 
apply in one and two dimensions where even the metallic 
state g > g c is poorly defined due to divergent Althshuler- 
Aronov and localization corrections to conductance at 
T = 0. The metal to insulator transiton/crossover at 
higher temperatures, T > gS in low dimensional metallic 
arrays was considered in Refs. [29ll3Qll3lj | . 

The temperature dependence of the conductivity of 
granular arrays which are insulators at T = (i.e. with 
g < g c ), is given by the activation law-^ 

a~a e- A ^ T , T < A M , (5) 

with <jo being the high temperature conductivity as long 
as the temperature is less than the Mott gap. Indeed, 
the finite temperature conductivity is due to the elec- 
trons and holes that are present in the system as real 
excitations . Their density is given by the Gibbs distri- 
bution that results in the conductivity dependence ©. 
Note that the expression JSJ) contains the renormalized 
Mott gap Am given by Eqs. (|2I3|) . 

It is important to note that at finite temperatures the 
continuous transition between the insulating and metal- 
lic phases, is, in fact, rather a crossover than the true 
phase transition. Indeed, the conductivity at nonzero 
temperature is finite in both phases, besides the Mott 
transition does not break any symmetry of the system. 
However, if the transition is discontinuous at T = it 
has to remain discontinuous in some finite temperature 
interval. In this case the Mott transition would be ob- 
servable even at nonzero temperature via the jumps of 
the physical quantities such as conductivity. 

The results presented above hold only for periodic ar- 
rays. Most of realistic granular /quantum dots arrays con- 
tain disorder, in particular, electrostatic disorder. 



B. Granular arrays with electrostatic disorder 

Electrostatic disorder causes fluctuations in electro- 
static energy of granules and can thus lift the Coulomb 
blockade at some sites of a granular sample. This re- 
sults, in its turn, in the finite density of states at the 
Fermi level and in appearing the variable range hopping 
as the dominating mechanism of conductivity. The bare 
density of states induced by the random potential can 
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be substantially suppressed due to the presence of the 
long-rang Coulomb interaction as it takes place in semi- 
conductors where the density of states (DOS) is given by 
the Efros-Shklovskii expressioniiSl 

v g {e) ~ (k/e 2 ) d \e\ d -\ (6) 

with e being the electron charge and R being the dielectric 
constant. This result was recently confirmed in Ref»2i 
(see also Ref^i) within the analytic approach of the lo- 
cator approximation. 

In the limit of weak coupling between the grains, g <C 
1, the model of granular array can be described in terms 
of the classical model that deals with the total charges 
on the grains only 

H ^fiifii + ^rii EfjUj, (7) 

i ij 

where \±i is the random external potential, is the elec- 
tron density on the site z, and Efj is the Coulomb in- 
teraction between the sites i and j. Taking into account 
the asymptotic behavior of the Coulomb interaction ma- 
trix Efj ~ e 2 /rijk, where r^- is the distance between the 
sites i and j and R is the effective dielectric constant in 
a granular sample we see that the classical model J7j| is 
essentially equivalent to the one studied by Efros and 
Shklovskii. Therefore their results must be applicable to 
the model of a granular array (J7J) as well. 

At the same time one expects that the DOS in an array 
of metallic granules is larger than that in a semiconductor 
since each metallic grain has a dense electron spectrum. 
Indeed, one has to remember that there are many elec- 
tron states in a grain that correspond to the same grain 
charge while in the model of impurity levels in semicon- 
ductors the charge is uniquely (up to the spin) identified 
with the electron state. The energy of the unoccupied 
state Si in the model (7J) is by definition the energy of an 
electron placed on this state. In the granular metal any 
state with energy larger than S{ will be also available for 
electron placement. Thus, in order to translate ES result 
(0 to the density of electron spectrum in granular metals 
one has to integrate the dependence (0 by the energy e 
and multiply it by the bare DOS in a single grain. As a 
result we obtain 

v(e)~vo(\e\R/e 2 ) d , (8) 

where v$ is the average DOS in a single grain (measured 
as the number of states per energy). However, the above 
expression cannot be used in the Mott argument for the 
hopping conductivity where one needs to estimate the 
distance to the first available state r within the energy 
shell e via the relation 




The problem with using the expression for DOS 0) 
in Eq. (0 is that expression (0 assumes strong corre- 
lations of the electron states in the space of the array 



coordinates. Indeed, expression 0) takes into account 
the fact that if there is a state available for placement 
of an electron with a given energy, then, typically on the 
same grain there will be plenty of other states available 
for electron placement. But, for application to the hop- 
ping conductivity we should not count different electron 
states that belong to the same grain since it is enough 
to find at least one to insure the transport. Thus, when 
finding DOS relevant for the hopping transport the low- 
est energy states within the each grain are to be count 
only. We arrive at the conclusion that while in granu- 
lar metals the electron DOS is modified according to 0), 
in order to find the distance to the first available state 
within a given energy shell via Eq.(0, even in granular 
metals, one has to use the ES expression for DOS in its 
form (0. Similar considerations were presented in Ref^. 
Following this reference we will call the DOS that counts 
only lowest excited states in each grain and is relevant for 
the hopping conductivity by the density of ground states 
and will denote by a subscript g as in Eq0 

Apart from the knowledge of the spectrum of electron 
excitations, description of the hopping conductivity re- 
quires understanding the tunneling mechanism of an elec- 
tron between the two distant grains with energies close to 
the Fermi level through an optimal chain of other grains 
where the Coulomb blockade is present. At low tempera- 
tures such tunneling process can be realized by means of 
the virtual electron hopping (cotunneling). As we men- 
tioned in Introduction one has to distinguish the elastic 
and inelastic cotunneling processes. The elastic cotunnel- 
ing is the dominant mechanism for hopping conductivity 
at temperatures T < ^/Eq5, while at larger temperatures 
the electron transport goes via inelastic cotunneling pro- 
cesses. 



1. Hopping via elastic cotunneling 

Considering the probability for an electron from the 
site zq to tunnel to the site in+i it is convenient to put 
the Coulomb interaction energy at these sites to zero and 
count the initial, £o> an d final, £jv+i? electron energies 
from the Fermi level. The presence of the electrostatic 
disorder on the grains is modelled by the random po- 
tential fjLi, such that the energy of the electron (hole) 
excitation on the site i is 

E±=E?±ih. (10) 

The probability P e \ of a tunneling process via the elas- 
tic cotunneling can be easily found for the case of the 
diagonal Coulomb interaction Efj = Ef Sij : 

P e i=wg N ^ (^-gj <S(6r+i-£o), (11) 

where the factor w = n(£o)(l — n(£/v+i)) takes into ac- 
count the occupations n(£o) and n(£/v+i) of the initial 
and final states respectively, and the bar denotes the 
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geometrical average of the physical quantity along the 
tunneling path. For example the average tunneling con- 
ductance g is defined as 



Ing 



1 



N 



N - 



i5> 



9k, 



(12) 



where the summation goes over the tunneling path and 
gk is the tunneling conductance between the k-th and 
(fc + l)-st grains. The energy E is the geometrical average 
lni^ = ^2k=i °f the following combination of 

electron and hole excitation energies 

E k = 2(l/E+ + l/E-y 1 . (13) 

The presence of the delta function in Eq. ((11)) reflects the 
fact that the tunneling process is elastic. 

The result (jTTJl for the tunneling probability P e \ implies 
that on average the probability falls of exponentially with 



the distance s along the path: P e i < 
dimensionless localization length ^ 



\n(Eir/cg5) 



e 2s /fez ? where the 
is given by 

(14) 



The numerical constant c in the above expression is equal 
to one for the model of diagonal Coulomb interaction, 



E tj = E ' 



Sij. We show below that the inclusion of the 



off-diagonal part of the Coulomb interaction results in 
the renormalization of the constant to some value 0.5 < 
c < 1. 

In a full analogy with the hopping conductivity in 
semiconductors^ in the case of granular metals inelas- 
tic processes are required to allow an electron to tunnel 
to a state with higher energy. In granular metals such 
processes occurs via phonons as well as via the inelastic 
collisions with other electrons; these proceses are not re- 
flected in the formula ((11)). Applying the conventional 
Mott-Efros-Shklovskii arguments, i.e. optimizing the full 
hopping probability oc exp[— (2s/£ e /) — (e 2 /RTas)\, one 
obtains Eq. (Q for the hopping conductivity with the 
characteristic temperature 



2 /ak^ e 



(15) 



where k is the effective dielectric constant of a granu- 
lar sample, a is the average grain size, and £ e z given by 
Eq. (when deriving JTSJ) we considered the tunneling 
path as nearly straight). 

In the presence of the strong electric field S , the direct 
application of the consideration of Ref. |32| for the case of 



low temperatures e , 
current dependence 



gives the nonlinear 



J ^ Joe 



-(£o/£) 1/2 



(16) 



where the characteristic electric field S is given by the 
expression 



el- 



(17) 



The results presented in this subsection are valid as long 
as the contribution of the inelastic cotunneling to the 
hopping conductivity can be neglected. This is the case 
at low temperatures and electric fields T, Sea <C ^/SE^. 
Below we present our results for the opposite case where 
inelastic processes give the main contribution to hopping 
conductivity. 



2. Hopping via inelastic cotunneling 

Inelastic cotunneling is the process where the charge 
is transferred by means of different electrons on each el- 
ementary hop. During such process the energy brought 
to a grain by incoming electron differs from the energy 
taken by the outcoming one. 

As in the case of elastic cotunneling we assume that 
the electron tunnels from the grain zq with energy £o to 
the grain ijv+i with energy £jv+i- For the probability of 
such a tunneling process through a chain of grains in the 
tunneling path we find the general expression 2 ^ 



4ttT TT^+i 



4ttT 



2N 



|r(jv + j^)| 2 * 

± n ± p 2T 

T(2N) 



(18) 



where T(x) is the Gamma function and A = £at+i — £o 
is the energy difference between final and initial states. 
The appearance of the factor e -A / 2T is consistent with 
the detailed balance principle. Indeed, at finite tempera- 
tures an electron can tunnel in both ways with either in- 
crease or decrease of its energy. According to the detailed 
balance principle the ratio of such probabilities must be 
e A / T , that is indeed the case for the function P in (HSJ), 
since apart from the factor e - A / 2T the rest of the equa- 
tion is even in A. 

To obtain the expression for the hopping conductivity 
one has to optimize Eq. (|T8j) with respect to the hopping 
distance N under the constraint 



NakA/e 2 - 1, 



(19) 



following from the ES expression for the density of ground 
states (0 . Optimization of Eq. ([T5|) is a bit more involved 
procedure than the standard derivation of Mott-Efros- 
Shklovskii law based on the Gibbs energy distribution 
function (see Appendix), however it leads to the essen- 
tially the same result, i.e. the ES law with 



T (T)~e 2 /ak^n(T), 



(20) 



with the dimensionless localization length ^ n (T) being 
weakly temperature dependent 



\n[E 2 /mcTrT 2 gY 



(21) 



where the coefficient c equals one for a model with diag- 
onal Coulomb interaction matrix. As in the case of the 
inelastic cotunneling, the long range interaction result in 
the reduction of the coefficient to some 1/4 < c < 1. 
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At zero temperature Eq. ([T5|) for inelastic probability 
Pi n for A < is simplified to 



Pzn(T = 0) 



w2 2N n I A 



2N-1 



(27V -1)! E 2N 



N+l 



) 



(22) 



while for A > Eq. ([T5|) gives zero since the process 
with increase of the energy of the tunneling electron is 
prohibited at zero temperature. 

To obtain the hopping conductivity in the regime of 
strong electric field £ and low temperatures, following 
Ref. 32], we use the condition (|TU|) that defines the dis- 
tance to the first available electron site within the energy 
shell A together with equation 



e£r ~ A, 



(23) 



that relates A with the electric filed £ and the distance 
between the initial and final tunneling sites r ~ Na. 
Together Eqs. and (J2HJ) define A ~ \/£e 3 /k and 
N ~ y/ e/k£a 2 that being inserted into Eq. ((221) res ult in 
the expression for the current 



J ^ Joe 



-(W£) 1/2 



(24a) 



where the characteristic electric field £q is a weak func- 
tion of the applied filed £ 



£o(£) 



Yn 2 [E 2 /e 2 £ 2 a 2 g}. 



(24b) 



Coulomb interaction matrix is related to the capacitance 
matrix in the usual way 



Ef. 



(e 2 /2) 07: 



(26) 



As we discussed in Introduction the random potential 
jii that models the electrostatic disorder removes the 
Coulomb blockade from some grains. 

We let the tunneling elements be random Gaussian 
variables defined by their correlators: 



(27) 



where i,j are the coordinates of the nearest neigh- 
bor grains. The intergranular conductance between the 
grains i and j is related to average matrix elements as 



(28) 



where v\ is the density of states in the grain i. The 
conductance gij is defined per one spin component, such 
that, for example, the high temperature ( Drude ) con- 
ductivity of a periodic granular sample is a = 2e 2 ga?~ d , 
with a being the size of the grain. In the following we will 
consider both cases of periodic ( in particular, it means 
that Hi = 0, Vi = vq, = g) and irregular arrays be- 
ginning with the former case. 



IV. MOTT GAP AT SMALL TUNNELING 
CONDUCTANCES 



III. THE MODEL 

Above presented results were derived within the follow- 
ing model: We consider a d-dimensional array of metallic 
grains. The electron motion within each grain is diffusive 
and the grains are assumed to be weakly coupled such 
that the intergranular tunneling conductance g is much 
less than the intragranule one. The system is described 
by the Hamiltonian 

H = H + H C + J2 t% kj #(*0 Mkj), (25a) 

ij 

where t% k , is the matrix element representing the elec- 
tron tunneling amplitude from the state ki of the grain 
i to the state kj of the neighboring grain j. The Hamil- 
tonian i^o in the r. h. s. of Eq. (|25a|) describes noninter- 
acting isolated disordered grains. The term H c includes 
the electron Coulomb interaction and the local external 
electrostatic potential fii on each grain 

H c = ^2,ni Efj fij + ( 25b ) 

ij i 



In the limit of small tunnelling conductances, g <C 1, 
and zero temperature a periodic granular array is a charge 
insulator. The key characteristic of the insulating state 
is the Mott gap Am in the electron excitation spectrum 
that can be defined as 

Am = £?a/=i — fi — E N =o, (29) 

where \i is the chemical potential, Ejy=o is the ground 
state energy of the charge neutral array and En=i is the 
minimal energy of the the system with an extra electron 
added to the neutral state. In the limit of vanishing tun- 
neling conductance, g — > 0, the Mott gap Am coincides 
with the Coulomb energy Eq = Ef i in a single grain 

Am g^o = (30) 

At finite intergranular coupling the Mott gap is re- 
duced as a result of the virtual presence of an electron 
on the neighboring grains. In the case of weak coupling, 
^ « 1, the resulting suppression of the Mott gap can 
be easily calculated within the perturbation theory. The 
first non vanishing correction appears only in the second 
order in tunneling matrix elements and can be found with 
the help of the textbook formula^ 



where Efj is the matrix of Coulomb interaction and hi is 
the operator of electron number in the i— th grain. The 



(31) 
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FIG. 2: This figure illustrates the virtual electron tunneling 
processes between neighboring grains that contribute to the 
reduction of the Mott gap. The tunneling processes shown on 
figure a) determine the correction to the ground state energy 
of the charge neutral state En=o- The correction to the energy 
of the ground state with one extra electron (placed on grain 
1) En=i is illustrated in figure b). 



where AE is the energy correction to the ground state 
and the matrix elements of the perturbation V are taken 
between the ground state and excited states k. 

Correction due to finite intergranular coupling should 
be included in both terms En=i and En=o in Eq. 
The ground state with one extra electron N = 1 in the 
limit of zero coupling is degenerate since the extra elec- 
tron can be present on any grain. Further, we will assume 
that it occupies grain 1 as it shown in Fig. 0b. Since we 
want to find the difference between the corrections to 
N = 1 and N = ground states, it is clear that within 
the second order perturbation theory in tunneling ele- 
ments tkxM on ly electron hops between grain 1 and its 
neighbors result in nonzero contribution. All other pos- 
sible hops give rise to equal corrections to the energies 
En=i and En=o that are mutually cancelled in Ea. (|^j) . 
Further, since all the neighbors of grain 1 are equivalent, 
we will consider only electron hops between grain 1 and 
one of its neighbors that we denote as grain 2. Con- 
tribution of the hops between grain 1 and all its other 
neighbors will simply result in the factor z (coordination 
number ) in the final answer. 

Now let us consider the correction to the ground state 
energy £?a/=o- In this case the matrix elements Vk,o cor- 
respond to the processes of electron tunneling between 
the neutral neighboring grains. Consider the tunneling 
process from grain 1 to grain 2 shown in Fig. The 
excitation energy of this process is Ek 1 +£k 2 + Eeh, where 
Sk 2 is the bare (with no Coulomb energy included) en- 
ergy of an electron excitation in grain 2, is the bare 
energy of a hole excitation in grain 1, and E e h is the 
Coulomb part of the energy of the electron-hole excita- 
tion E e h = E\ x + E 22 — 2£^2 that for the periodic case 
under consideration is reduced to 



eh 



2EZ - 2E 



12- 



(32) 



The energy correction corresponding to such process is 

I^!,/C 2 | 2 



-AE N=0 = 2 

/Ci,/C2 



(33) 



where the factor 2 takes into account the similar process 
of electron hopping from grain 2 to grain 1. 

Analogously we find the correction to the energy Ejy = i 
as illustrated on Fig. 03: The excitation energy of the 
process of electron tunneling from grain 2 to grain 1 is 
Ski + e k2 + 3JS7f! + E c 22 - 4Ef 2 = e kl + e k2 + 2E eh , while 
the excitation energy of the process of electron tunnelling 
from grain 1 to grain 2 is Sk x +£k 2 ^E\ x —E 22 = £ k ± +£/c 2 - 
The corresponding correction to the energy En=i reads 



-AE 



N=l 



ki ,k 2 



\t kl , 



2E, 



eh 



£fci + £k 2 



(34) 



One can see that corrections (|33l l3l|) are ultraviolet di- 
vergent. However, their difference must be finite since 
it represents the measurable quantity - correction to the 
Mott gap. It is, indeed, the case: Subtracting Eq. 
from Eq.JHIJ), presenting summation over states as inte- 
grals £fci(fc 2 ) —> e i(2) an d introducing further the variable 
e = E\ + £2 we obtain 



2.2 

ZVnt 



ede 



1 



2E ( 



eh 



where is the density of states in a single grain and t 2 
defined by Eq.JSTJ) appears from the averaging over ma- 
trix elements. Taking the integral in the above expression 
we obtain the result (J2J). 

In the derivation of the reduction of the Mott gap pre- 
sented in this section we neglected the fact that an extra 
electron added to the neutral system in the presence of fi- 
nite intergranular coupling can move diffusively over the 
sample. Contributions corresponding to such processes, 
as we will see below, are suppressed by an extra small 
factor S/Eq <C 1 and, thus, can be neglected. 



V. PHASE AND CHARGE REPRESENTATIONS 

The perturbation theory presented in Sec.|^is limited 
to the regime of low tunneling conductance g <C 1. More- 
over, the high order tunneling processes that are very im- 
portant for understanding the physics of Mott transition 
and hopping conductivity cannot be included in this way. 
For this reason one has to seek for a convenient approach 
applicable at intermediate tunneling conductances and 
capable of description of the effects of high order tun- 
nelling processes. In this view, a useful tool in approach- 
ing this problem is the phase functional approach^ that 
allows further mapping of the system to the Coulomb gas 
of charged particles in d+1 dimensionsiiSi. 
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A. Phase representation 

The phase functional approach is based on the de- 
coupling of the Coulomb interaction term ([25b[) in the 
Hamiltonian (|25a|) by means of the the potential field 
V(r) that is further absorbed into the gauge transforma- 
tion of the fermionic fields 

Mr) ^Mr)e~ iMr \ &(r) = f V^dr' . (35) 

J — oo 

Since the action of an isolated grain is gauge invariant 
the phases <pi(r) enter the whole lagrangian of the system 
only through the tunnelling matrix elements 

tij > tij e ^ 1 ^ \ (36) 

where 4>ij(r) = 4>i{r) — (j)j(r) is the phase difference of 
the i— th and j— th grains. 

A simple description in terms of phase variables can 
be obtained via integrating out the fermionic fields and 
expanding the resulting rather complicated action to the 
lowest non vanishing order in tunneling elements. As a 
result one obtains the Ambegaokar, Eckern, Schon (AES) 
actioni& 



Sa 



ES 



-y 

2e 2 ^ 



dr I 



• ijii) Cij (<j>j + ijij) + S t [4>], 



(37a) 



where the tunneling term St at zero temperature is 



= 2^ E 9ij 



<IJ> 



a i(f)ij (ri)—i(f)ij(T2) 

(ti-t 2 ) 2 



(37b) 



where the angular brackets mean summation over the 
nearest neighbors. The higher order terms in the ex- 
pansion of the total action in the tunneling matrix ele- 
ments cannot be always neglected and, in general, must 
be taken into account along with the nearest neighbor ac- 
tion (|37|) . However, before considering these more com- 
plicated higher order terms, it is instructive to analyze, 
first, the action ([57| in details. Moreover, as we will 
show, even at T = there is a regime where the AES ac- 
tion is applicable. We, thus, proceed to the description 
of the charge mapping beginning with the AES action, 
and later, in Sec. lVIlfl we will show how the higher order 
tunneling terms can be included into the mapping. Also, 
for simplicity, we will first consider the case of the regular 
periodic array. 



B. Charge representation 

Proceeding to the mapping of the action ([57)1 onto the 
model of the classical Coulomb gas and using the method 
of Ref . ^3 •> we fi rs t expand the partition function Z to 
all orders in the tunneling part of the action St 



■/AH. 



(38) 



Here Zq is the partition function of a system of isolated 
grains and the angle brackets assume averaging with re- 
spect to the Coulomb phase action that for the case 
\±i = under consideration reads 



S c 



—y 

2e 2 ^ 



dr d>i d 



(39) 



On the next step, for each term in the expansion (|38|) we 
write the N-th power of the tunnelling action St explic- 
itly as the product 



N 



(n) 



(40a) 



where the index n of the tunneling action St reminds that 
the time integrations in each term in the above product 
are independent 



4 ~ 27T ki (r -- r « ; 



a i( t>ij(.r n )-i(f>ij(r^) 



(40b) 



Writing the powers of the tunneling action as products 
([40 a)) allows to implement the integrations over the phase 
fields (j)i{r) exactly for each term in the expansion ([55)l 
with the help of the formula 



where U c given by the expression 



n\ 7~ n <2 | 6 ni 6 n2 i 



(41) 



(42) 



can be viewed as a Coulomb energy of a system of 
classical charges in a space with an extra time dimen- 
sion^. The classical charges e n = ±1 interact via the 
Id Coulomb potential along the time direction while the 
strength of the interaction is given by the "real" Coulomb 
energy of original quantum electrons E\- . 

From Eq. ([40b[) we see that each tunneling term S^ [4>] 
contains the sum of four phase fields in the exponent 



iM r n) ~ ^i( r n) - i<t>i{r' n ) + i(j>j{r' n ). 



(43) 



Each of these fields, according to Eqs. (|41I42|) . gives rise 
to a classical charge located at coordinates corresponding 
to the arguments of the phase fields. This means that the 
classical charges appear as quadruples that occupy only 
the neighboring sites as it is shown in Fig. [51 

The time dependent kernel in Ea. ([40b[) can be written 
in the form 



2tt (r n - r' n f 



-WjT n -T> n ) 



(44) 



such that defines the "internal" energy of n— th 
quadruple 



AT=1 



^ = ln[27r(Ar„) 2 /ff], 



(45) 
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FIG. 3: Illustration of the mapping of a quantum granular 
array model in the nearest neighbor hopping approximation 
onto the classical model of gas of quadruples in d + 1 dimen- 
sions. 

with Ar n being the size of the n— th quadruple in the r— 
direction as shown in Fig. 03 

Thus, the partition function, Eq. (|38|) . of the quan- 
tum problem under consideration in the nearest neigh- 
bor electron tunneling approximation can be presented 
as a partition function of classical charges that appear 
in quadruples all charges being subject to the Coulomb 
interaction ((42)). The total energy of the system of N 
quadruples can be written as 

N 

U* = Y,K + Ut N , (46) 

n=l 

where U% is the internal energy of the n— th quadruple 
and U% N is the Coulomb energy of 4N charges that form 
N quadruples defined by Eq. P^|) . We notice that dif- 
ferent quadruples interact with each other only via the 
Coulomb part of the energy U c . 

In the following sections we will consider the thermo- 
dynamic properties of the classical model and set corre- 
spondence between the parameters of the original quan- 
tum model and the effective classical one. 



VI. THERMODYNAMICS OF THE CLASSICAL 
MODEL IN THE LIMIT OF LOW QUADRUPLE 
DENSITY 

The thermodynamic potential Q of the classical 
model can be found in the limit of the low quadru- 
ple density that, as we will show below, corresponds to 
the case of the low tunneling conductance g <C 1 in the 
quantum model. In this limit one can neglect the mu- 
tual interaction of deferent quadruples considering them 
as free independent objects and the thermodynamic po- 
tential Q of the whole system is reduced to the sum of 
the free energy JF of a single quadruple multiplied by 
the number of the quadruples N and the entropy S of 



the quadruple configuration space 

n = NFo - S. (47) 

The entropy of N independent quadruples is given by the 
expression 

S = N\nV -\nN\, (48) 

where V = LM is the volume of the d + 1 dimensional 
system with L being the size of the system in r— direc- 
tion and M being the total number of the sites (grains) 
in the array. The term In TV! in Eq. ([48)) appears from 
the term 1/N\ in the expansion ((38)). In the classical sta- 
tistical mechanics such term appears due to equivalence 
of the classical configurations that differ only by mutual 
permutations of particles. Thus, in the thermodynamic 
limit N —> oo for the thermodynamic potential Q in the 
case of low quadruple density we obtain the expression 

n = NJ 7 o + Nhip-N, (49) 

where p = N/V is the density of quadruples in d + 1 
dimensions. 

The partition function of a single quadruple is given 
by the expression 

z«=- r% e ~ EehT i ( 5 °) 

where the energy of the electron-hole excitation is denned 
by Eq. (|52|) . Equation (|5U|) is divergent on low integration 
limit and in order to regularize it we introduce the time 
cut-of A obtaining Zq = g/irX so that the free energy of 
a single quadruple T§ — — In Zq becomes 

^o = -M0ArA). (51) 

A. Density of classical charges 

The thermodynamic potential given by Eas. ()49l51)) de- 
termines all the thermodynamic quantities of the classical 
model in the limit of low quadruple density. In particu- 
lar, the density itself can be obtained via minimization of 
ft with respect to the number of quadruples TV resulting 
in 

P = 0ArA. (52) 

We would like to note that the obtained density coincides 
with the single quadruple partition function p = Zq. The 
result for the density p (|52|) can be better understood 
considering a discrete time axis. In this case the cut-of 
distance A would be the distance between the neighboring 
sites in r— direction and the dimensionless density ac- 
cording to Eq. (|52|) would be g/ir. Thus, it becomes clear 
that the limit of low conductance g <C 1 corresponds 
to the limit of the low occupation number of classical 
charges. 
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B. Relation between the Mott gap and dielectric 
constant of classical gas 

In order to make the developed mapping useful one has 
to express the physical observables of the original quan- 
tum model in terms of the corresponding thermodynamic 
quantities of the classical model. Below we show that the 
key characteristic of the insulating state - the Mott gap 
Am is related to the dielectric constant n T of the effec- 
tive classical model measured in the direction of r— axis 
as 

A M = E c Jn T . (53) 

To prove relation (|53j) we first note that the gap Am 
in the electron spectrum determines the asymptotic be- 
havior ( r — > oo ) of the single electron Green function 

Gu(t) ~ e~^ M ' T '. (54) 

At the same time the asymptotic behavior of the Green 
function Gu (r) with exponential accuracy coincides with 
the asymptotic behavior of the phase correlation function 

Gii{ Tl - r 2 ) ~ $(n - r 2 ) = ( e ^(^)-^(^) ), (55) 

where the angular brackets stand for averaging with re- 
spect to the phase functional (|37[) . written explicitly 

*(ri - r 2 ) = Z' 1 J D[<t>] e^^)-^(^) e S[<^ ^ 

All the considerations that was previously used in order 
to derive the effective model can be applied to the 
correlation function (|56|) as well. Repeating all the nec- 
essary steps we notice that phases <pi (n ) and & (t 2 ) that 
enter the exponent in Ea.(|56j) result in appearance of 
two external charges of opposite signs located on the site 
i at time coordinates t\ and r 2 . Further, we see that the 
phase correlation function <£> is related to the thermody- 
namic potential of the system with two external charges 
^2<Vi - r 2 ) as 

*(r) =e-^ 2(r) . (57) 

The function £l 2 (r) is normalized such that ^ 2 (0) = 0. 
The function Sl 2 is nothing but effective interaction of 
two external charges in the presence of quadruples. The 
initial electrostatic interaction is reduced in the presence 
of quadruples and the asymptotic behavior of the func- 
tion Q2 takes the form 

« 2 (t) = M^o/«t, r^oo, (58) 

where n T is the dielectric constant. This proves the rela- 
tion between the Mott gap and dielectric constant of the 
classical Coulomb gas in r— direction ((551 ). 



a) 



b) 
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FIG. 4: An electric filed that charge e<z feels consists of the 
field created by the charge e\ and the fields created by the 
charges of the quadruple qi, ...q±. 

C. Mott gap in the limit of gz <C 1 

Finally, in order to better understand the relationship 
between the original quantum model and the effective 
classical one we will rederive the result (0) for the Mott 
gap in the case of small tunneling conductance « 1 
within the framework of the effective classical model. 

Let us consider two charges e\ and e 2 ( see Fig. 0| ) 
of the opposite signs located on the site i at time co- 
ordinates Ti and r 2 respectively. In the absence of sur- 
rounding quadruples their interaction energy would be 
Eq\ti — r 2 | that gives the zero-order value of the Mott 
gap Am = Eq. The bare interaction energy Eq of two 
charges is reduced due to the presence of quadruples. In 
order to find this reduction it is convenient to consider 
the electric field created by the charge 1 and partially 
screened by the surrounding quadruples at the point r 2 
of location of charge 2. In the limit of low quadruple 
density the screening of the interaction due to different 
quadruples can be considered independently. 

It is easy to see that the electric field at the point r 2 
is affected by the presence of the quadruple only if the 
later is placed such that charge e 2 is inside the quadruple 
as it is shown in Fig. 2] Depending on the arrangements 
of quadruple charges one can distinguish two different 
situations shown in Fig. 2] a) and b). In the case a) the 
electric field that charge e 2 feels is 

E^ = -3£ C + 2E C 12 = -E C Q - E eh , (59a) 

while in the case b) it is 

= E c - 2E C 12 = -E c + E eh . (59b) 

We see from Eqs. ((55)) that the electric field at the point 
r 2 of location of charge 2 gets an additional contribu- 
tion ±E e h depending on the arrangements of quadruple 
charges. The probability of the configuration a) is given 
by the Gibbs weight as 

where the time coordinates of quadruple charges r qi and 
T q2 are taken with respect to coordinate r 2 where charge 
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e2 is located such that r q2 > and r qi < 0. Analogously 
we find the probability of the configuration b) 

Pb (r qi ,T q2 )= 1 ^ -^1 (^2 ~r qi ) ~r qi E eh 

The probabilities P a and must be further multiplied 
by the quadruple density p, coordination number z and 
combinatoric factor 2. Multiplication by p simply can- 
cels the factor Z^ 1 due to the relation p = Z$ that we 
pointed out below Eq.(|52J). The additional contribution 
to the electric field at coordinate T2 of charge 2 due to 
the presence of the quadruples becomes 

/0 poo 
dr ql / dr q2 (P b - P a ). (60) 
-oo JO 

The dielectric constant n T is given by the expression 
n' 1 = (Eft - AE)/Efi. Taking integrals in Eq. {EJ we 
finally get the dielectric constant 

«r =1 ^T ln2 ' ( 61 ) 

7T £/q 

that taking into account relation ((55)) between the Mott 
gap and dielectric constant is indeed consistent with re- 
sult (0) previously obtained within the direct perturba- 
tion theory in tunneling matrix elements. 



VII. MOTT GAP AT LARGE TUNNELING 
CONDUCTANCES 

In this section we continue our investigation of the 
quadruple model ((IS)) now coming to the regime of mod- 
erately large tunneling conductance gz > 1. We show 
that within this model the Mott gap Am is exponentially 
reduced with growth of g but nevertheless it always re- 
mains finite and, thus, the quantum system represented 
by the model ((IS)) is always an insulator at T = even 
at arbitrary large tunneling conductance g. This state- 
ment, however, turns out to be an artifact of neglecting 
the high order tunnelling terms in derivation of the phase 
action ((57)) and the following from it model ((IS)) . Below 
we show that models (|37I46|) are applicable as long as the 
Mott gap Am is larger than the inverse escape time of 
an electron from a grain T = gS. Thus, the results of this 
section are applicable as long as the obtained Mott gap 
Am is larger than the inverse escape time g5. 



A. Mean field approximation 

The model ((IS)) can hardly be solved in general case 
and further in this section we make two assumptions 
that will allow us to construct a mean field approach^. 
Namely, 

(i) we take the diagonal matrix of Coulomb interaction 
Ef j =ESS ij . 



Ql® © Q2 






V(t) 



FIG. 5: Screening of the external charge e in the dipole 
model. The Coulomb potential of the charge e is screened by 
the dipole charges q\ and q^ resulting in the effective screened 
potential V(r). 

(ii) we assume that the coordination number z is large. 
These two assumptions allow to reduce the quadruple 
problem in d + 1 dimensions to the problem of one di- 
mensional dipoles with self consistent dipole interaction. 

To construct a mean filed theory let us pick up a site 
i and consider the influence of the neighboring sites on 
it in a mean field manner. The quadruple that occu- 
pies the site i and a neighboring site j is reduced to the 
dipole located on the site i. The mutual interaction of two 
charges of the resulting dipole consists of the quadruple 
interaction U q and of the Coulomb interaction energy of 
quadruple charges located at the site j. The later must be 
taken into account in a self-consistent way. The mutual 
influence of the neighboring sites due to direct Coulomb 
interaction is absent due to our choice of the diagonal 
Coulomb matrix Efj. Thus the internal dipole interac- 
tion energy becomes 

U d = \n[n(Ar n ) 2 /gz] + tt 2 (Ar n ), (62) 

where the interaction energy f^fV) of two external 
charges must be calculated self consistently in terms of 
the dipole model. The asymptotic behavior of this func- 
tion at the same time according to Eq. ((55)) determines 
the dielectric constant n T that is directly related to the 
Mott gap Am of the original model. 

Thus we arrive at the simplified one dimensional model 
of dipoles bounded by the internal interaction ((52)). 
Apart from the bounding interaction all charges are sub- 
ject to the 1-d Coulomb interaction 

U id = -y l Tn i r n 2 | e n2 . (63) 

ni ,ri2 

In the limit of large density p ^> 1 that corresponds 
to the case of large tunneling conductance g ^> 1 the 
Coulomb interaction gets sufficiently screened. For this 
reason the expression for the density ((52)) as well as the 
identity Zq = po are applicable in the dense limit too. In 
the dense limit the screening effects can be usually well 
described within the Debay theory: 

Consider the charge e placed at r = shown in Fig. 
The Coulomb potential of this charge screened by the 
dipole charges obeys the Poisson equation 

d 2 T V(r) = -2E< (q(r) - 8(r)), (64) 
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where q(r) is the average density of all dipole charges and 
S(t) is the delta function. At the same time the density 
of dipole charges is given by the Gibbs distribution: The 
contribution of two charges q\ and q 2 of a single dipole 
( located at coordinates t\ and r 2 respectively ) to the 
charge density q(r) is 

^o" 1 [S(r - n) - S(r - r 2 )] e -^(-i-^)-^(n)+V(r 2 ) ? (65) 

where the energy Ud(j) that enters the Gibbs expo- 
nent includes the dipole interaction energy Ud(r) and the 
screened Coulomb interaction of charges qi and q 2 

U d (T)=U d (T)+n 2 ( T ). (66) 

We note that within the approximation that we use 
the screened potential V(r) is nothing but the effective 
screened interaction of two charges Vt 2 (t) so that 

V(r) = n 2 (r). (67) 

Multiplying Eq. (|65[) by the dipole density p, making use 
of the identity Zq = p and integrating over the coordi- 
nates T\ , T2 we obtain the averaged charge density 

2az f e -2V(T-r') 

Q(r) = J - / dr'- — sinner') - V(r)}. (68) 

7T J (r - t'Y 

— 00 

Equations (|64[) and (|68|) form a closed system of equa- 
tions that determine the potential V(r). In the next sub- 
section we present the solution of these equations. 



B. Solution of the mean field equations 

Numerical solution of Eqs. (|64I68|) for moderately large 
tunneling conductances gz > 1 shows that the external 
charge is never screened completely and the long time 
asymptotic behavior of the function V(r) in Eq. (|67|) has 
a linear form 

V(t) = \t\E c q k-\ t^oo, (69) 

where the inverse dielectric constant k~ x becomes very 
small for large conductances but nevertheless it remains 
finite. The numerical solution also shows that the poten- 
tial V(r) behaves logarithmically on intermediate time 
scales and then turns to the insulating linear behav- 
ior on longer times (crossover scale is given below). 
This suggests the following analytic approach to the so- 
lution of mean field Eqs. (|64I68|) : To find the potential 
V(r) on intermediate times where nonlinear effects are 
not important we consider Eq. (|68|) in the linear approx- 
imation with respect to V(r). As a result in the Fourier 
representation we obtain 

q u = -2zg\uj\V u . (70) 



Together with the Poisson equation written in the 
Fourier representation — uo 2 V(uo) = —2Eq (q u — 1) equa- 
tion (f70j) defines the potential 

V " = ~u) 2 /2E% + 2zg\u)V 

that being transformed to the time representation for r > 
(E^g) -1 results in 

V(r) = (72) 

The linear approximation that we used in Eq. is 
valid as long as the potential V(r) is small. This assumes 
that dependence (|72|) is applicable for not very long times 
r <C exp(7rzg)/ gE$. From the logarithmic form of the 
answer (J72|) it follows that the change in the charging 
energy 

E c ^ZE c (73) 

in Eq. (|T2|) results in a constant shift of the potential 
V(r) — > V(r) + (\j2i\zg) ln£. At the same time the only 
term in Eq. (|68j) that is sensitive to this shift is the po- 
tential V(t) in the exponent that provides the renormal- 
ization of conductance 

5 - 5 -(l/7rz)ln£. (74) 

This consideration allows to solve the system of 
Eqs. (|64I68|) for the case of large tunnelling conductance 
in the renormalization group(RG) scheme. Since the 
charging energy Eq is the only energy scale in the prob- 
lem, the scaling of this energy with conductance deter- 
mines the scaling of the Mott gap. Thus, from the scaling 
Eqs. (| 7 31 74)) we finally obtain the expression for the Mott 
gap ©. 

Equation (0 agrees with the result of RG approach 
of Ref. Q where the opening of the Mott gap Am was 
assumed to happen at the point where the tunneling con- 
ductance running under RG scheme flowed to values of 
the order one. We would like to note that the present 
approach does not rely on any assumptions of this kind: 
we first apply the RG scheme to reduce the conductance 
to values gz > 1, and further we solve the equation nu- 
merically and show that opening of the gap indeed takes 
the place, i.e. the asymptotics of the potential V(r) has 
a linear form indeed. 



VIII. INCLUSION OF HIGHER ORDER 
TUNNELLING TERMS 

So far we considered the quadruple model (fl?)|) that 
includes only the nearest neighbor electron tunnellings. 
In the present section we describe how the high order 
tunneling terms can be included into our mapping scheme 
and in the following sections we apply this technique to 
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FIG. 6: Diagram a) represents the lowest order correction in 
tunneling conductance g to the partition function, the first 
term in the expansion ((38)) . The solid lines denote the prop- 
agator of electrons, the crossed circles stand for the phase 



dependent tunnelling matrix elements Uje 



Diagram b) 



describes a higher order correction in tunnelling matrix el- 
ements Uj. The dotted lines represent impurity scattering. 
The diagram c) shows the zero dimensional diffusion in time 
representation. 



the problem of hopping conductivity and Mott transition 
in granular arrays. 

Keeping in mind the applications to the irregular ar- 
rays, as in the case of hopping conductivity, in this sec- 
tion we include the random on-site potential \i{ in the 
Hamiltonian into consideration. Also we take into ac- 
count variations of tunneling conductance g from grain 
to grain as well as variations of charging energies Ef and 
the mean energy level spacings Si. 

In order to construct the effective classical model that 
takes into account higher order tunnelling terms we pro- 
ceed in the same way as in the derivation of AES func- 
tional Eqs. (pTf|). First we decouple the Coulomb inter- 
action term in Eq. ()25aj) with the help of potential field 
V{t) 



1 

2^2 



Y^Vi + ifjii] dj [Vj + ifij] - i^niVi, (75) 



where rii = rii(r) is the electron density field. The term 
—iriiVi can be absorbed by the fermion gauge transforma- 
tion given by Eqs. ((55)) such that the action that governs 
dynamics of the potential field V(r) becomes 



l 

2^2 



Y^iVi+wWijM+Wj]- ( 76 ) 



Further, following Ref. [36| we proceed with diagram- 
matic expansion of the partition function Z in tunnelling 
matrix elements Uj e l ^ ij ^ . In the lowest order we repro- 
duce the results obtained within the phase functional. In- 
deed, the lowest order correction to the partition function 
Z shown in Fig. Eh simply results in the first term in the 
expansion JSHJ). All disconnected diagrams of the type 
shown in Fig. EJi must also be included in the partition 
function, this reproduces all high order terms in the ex- 
pansion (|38j). Integration over the phase fields </>i(r) can 
be implemented exactly as it was described in Sec. IV Bl 
resulting in the appearance of classical charges in d + 1 



T 




FIG. 7: Illustration of the mapping of the model of quan- 
tum electrons in a granular array onto the model of classical 
charges imbedded into the electron world lines forming loops. 
The illustration is shown for the case of 1 + 1 dimensions (Id 
granular array) , generalization to higher d + 1 dimensions is 
straightforward. 



dimensions with the Coulomb interaction energy 
U= \ ^irj^i-^kn! e n2 +^/i in r n e n . (77) 

ni,ri2 n 

that differs from Eq. only by the presence of the 
random potential /^ not previously taken into account. 
Thus, inclusion of all disconnected diagrams, Fig. Eh, and 
their further averaging over the phases reproduces the 
classical quadruple model ((46)). The expression for the 
internal quadruple energy for the case where conductance 
g varies from contact to contact is the straightforward 
generalization of Eq. ((15)) 



W n = ln[27r(Ar„) 2 /5i. 



(78) 



The advantage of the mapping of the original quan- 
tum model on the classical electrostatic system is that 
it allows to include the higher order tunneling processes 
shown in Fig. UJ) in essentially the same way as the near- 
est neighbor hoppings. The higher order diagrams in- 
clude the single grain diffusion propagator 



D = 27T(T7|cji -u 2 \ 



(79) 



that being transformed into the time representation re- 
sults in the expression shown in Fig. HJ. Each tunneling 
matrix element in the diagram Fig. [(J) includes the phase 
variables as tije l ^ ij ^ T \ This allows to present the diagram 
in Fig. IHb as a charge loop shown in Fig. [7| The charges 
interact electrostatically and are subject to the local po- 
tentials according to Eq. ((77)). The internal energy of a 
single K— th order loop is 

K-l 



EW = \n[27rK/go}-J2Hg k S/27T}+ln\AT \ 

k=l 
K-l 

+ ln|ArK|+^ln(|Ar fc | + |A^|), (80) 



k=l 
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where the time intervals Ar^ and Ar' k shown in Fig. [7| 
are restricted to have opposite signs, i.e. Ar^Ar^ < 0. 

We mapped the model of quantum electrons in a metal- 
lic array with Coulomb interaction onto the classical 
model of Coulomb charges bounded by the electron world 
lines and forming loops. First order loops that form 
quadruples represent the process of virtual electron tun- 
nelings to the neighboring grains and back. Higher order 
loops represent electron motion on larger distances. 

In the Mott insulator state electrons cannot move on 
large distances. This assumes that in the classical repre- 
sentation of the Mott insulator state the thermodynamics 
of the classical charges will be dominated by the contri- 
bution of the short loops, while high order loops will be 
energetically suppressed. 

In the metallic state, on the contrary, diffusive motion 
of electrons on any distances is possible meaning that in 
the classical model the loops of infinite order are present. 
It is clear that the insulator to metal transition in terms 
of the classical model corresponds to the transition at 
which the infinite loops that were suppressed in the insu- 
lating state begin to proliferate. Infinite loops represent 
almost free charges that screen the Coulomb interaction. 
The inverse dielectric constant k~ x of this state is zero 
that corresponds to the absences of the Mott gap Am- 

The detailed description of the Mott transition is diffi- 
cult since it requires the analysis of contribution of high 
order loops in the regime where the tunneling conduc- 
tance g is not small and the Mott gap Am is already 
strongly reduced due to the quadruple screening accord- 
ing to Eq. (PJ). Such description, on the qualitative level 
is performed in Sec. EJ while in the next section we will 
study the high order electron tunnellings in a more sim- 
ple regime of low tunneling conductance. As we discussed 
above, the high order tunnelings being not important for 
the thermodynamics in the regime of small g are crucial 
for the description of the hopping transport in the case 
of irregular array. 



IX. HOPPING CONDUCTIVITY IN 
GRANULAR METALS 

In Sec. IVIIII it was shown that the periodic granular 
array with low tunneling conductances g is the Mott in- 
sulator with the hard gap Am in the electron excitation 
spectrum. This assumes the activation conductivity be- 
havior, a ~ exp (— Am/T), and excludes the possibility 
of hopping transport. However, as we pointed out in the 
Introduction, in realistic granular arrays there are many 
reasons for electrostatic disorder that compensates the 
Coulomb blockade on some sites. This gives rise to the 
finite density of states on the Fermi level and opens the 
possibility for the variable range hopping transport. 

In this section we discuss the hopping conductivity in 
granular metals via elastic and inelastic electron cotun- 
neling. We first concentrate on hopping via elastic co- 
tunneling and later discuss the inelastic mechanism. 




FIG. 8: The electron world line representing the probability 
of elastic cotunneling from the 0— th to the N + 1— st grain. 

A. Hopping via elastic cotunneling 

To analyze the hopping transport mechanism in gran- 
ular array one has to consider the probability of electron 
tunnelling between the two distant grains with energies 
close to the Fermi level. At very low temperatures such 
a tunneling can be realized only as a virtual tunneling of 
an electron through several grains being in the Coulomb 
blockade regime - such a process is called the elastic co- 
tunneling. 28 

In the following we consider the probability of elastic 
cotunnelling through several grains within the approach 
of effective classical model. We assume that the elec- 
tron tunneling takes place between two sites zq and in+i 
where the local Coulomb gap is absent. In our model 
it means that the Coulomb energy on these two sites is 
compensated by the external local potentials \±i. Thus, 
removing the charge from the grain zq as well as placing 
it on the grain zn+i does not cost any Coulomb energy. 
The probability of tunnelling between these two states 
can be found as follows: In the basis of the eigenstates 
of a system of isolated grains the amplitude of tunnelling 
process from the site zq to the site zn+i is 

A ioM+1 =<0\c iN+l Scl o \0>, (81) 

where S is the evolution operator that written in the 
interaction representation takes into account only tun- 
nelling part of the Hamiltonian (|25aj) . Consideration of 
the amplitude Ai 0j i N+1 , however, is not convenient since 
the amplitude becomes zero after averaging over tun- 
nelling matrix elements tfci,fc 2 - This happens because the 
phase of the amplitude of tunnelling process is very sen- 
sitive to the particular realization of disorder. For this 
reason one has to consider the probability of the tun- 
nelling process P = A* A. The conjugate amplitude A* 
can be written as the probability of the inverse process 
of tunnelling between zjv+i— st and zq— th grains 

A* , iN+1 = A iN+ul0 =< 0| c i0 S 4 +i |0 > . (82) 

The probability of the tunnelling process P = A*A can 
be presented as an electron world-line shown in Fig. |S1 
The first interval of the world line (t^to) represents the 
initial state of the system where the electron is located at 
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the site zq. The following part ( tq, tat+i ) represents the 
amplitude of the tunnelling process A 



^o,^N+l 



from the site 

io to the site zjv+i- The interval (tjv+i,t|^ +1 ) describes 
the system with the electron located on the site i jv+i ■ The 
next interval (rj v+1 , Tq) represents the amplitude Ai 7V+l5 i 
of the process of tunnelling back to the site io, and the 
last interval (tq,tj) represents the system back in the 
initial state where the electron is located at site zq. 

The world line that represents the tunnelling probabil- 
ity P = A* A is very similar to the closed loops in Fig.[7| 
The only difference is that in the world line of Fig. [SI the 
Green functions on the sites io and iN+i must be taken 
in a different way: In order to capture the energy reso- 
lution of electron energies of the initial and final states 
we will assume that both sites io and zjv+i have only one 
quantum state each with energies £o an d £at+i (counted 
with respect to the chemical potential) respectively. 

When electron tunnels from site zq to site zi, the hole is 
created at site zq. Such process is described by the Green 
function of a negative time argument 37 



Go{tq - t ) = n(£o) e«°^- T °\ 



(83) 



where n(£) is the electron occupation number that re- 
minds that the electron is assumed to be initially present 
at site zq. 

Similarly, the process when the electron is located on 
site iN+i is described by the Green function of positive 
arguments 

Gn+M+i-tn+i) = [l-n(^+i)] e-^^+i-^+O., 

(84) 

where the occupation number 1 — n(£/v+i) insures that 
the final state must be empty for an electron to be able 
to tunnel there. 

Before proceeding to the calculation of tunneling prob- 
ability let us first fix notations for all time coordinates in 
Fig. [SI To distinguish the time coordinates that belong 
to the left part of the world line representing the direct 
tunneling process from the ones that belong to the right 
part representing the inverse process we denote the later 
ones by primes. Further, to distinguish the time coordi- 
nates of positive and negative charges we use the index 
+ and — respectively. 

The tunneling probability P represented by the world- 
line shown in Fig.[S|can be calculated easily for the case of 
a diagonal (short range) Coulomb potential Ef- = SijEf. 



1. Short range Coulomb interaction 

In this case time integrations on each site k can be im- 
plemented independently and the tunneling probability 
P is given by the product of probabilities P k representing 
the single site contributions. To find a single site contri- 
bution P k we note that there are two possible arrange- 



ments of the tunneling times: (i) r k - < r k + < r' k _ < r } 
and (ii) r k+ < r k - 



< T U < T 'k- 



k+ 

Let us first consider 



the case (i): The Coulomb energy of such configuration 
is given by 



u. 



E+ (Ar k + At£), 



(85) 



where the electron excitation energy E k is given by 
Eq. ((TJJ and the time intervals Ar k , Ar k > are defined 
as Ar k = r k + - r k - , Ar k = r k+ - r k _ . The internal 
energy of the configuration according to Eq. ([80)1 is 

U k = - Hg k S/2ir} + ln(Ar fc + Arj£), (86) 

where g k is the tunnelling conductance between the k— th 
and (k + 1)— st grains. One can see that the second ar- 
rangement of charges (ii) differs only by the expression 
for the Coulomb energy 

^=^-(|Ar,| + |A^|). (87) 

Now with the help of Eqs. ESI EH) we find the el- 
ementary probability P k determined by the sum of the 
contributions (i) and (ii) as 



Z7T 



e k 



E+(At+At') + e -E-{Ar+Ar') 



At + At' 



and taking the integrals we obtain 

P _ 9k S 
h ~ — f - • 
irE k 



(88) 



where the energy E k is the combination of the energies 
E+ and E~ defined by Eq. f TS ) ). 

To find the time dependence of the tunneling process 
we notice that the time intervals Tq — To and t' n+1 — 
tat + i coincide within the accuracy of the inverse charging 
energy. In the limit that we consider the charging energy 
is the largest scale in the problem that allows to take both 
intervals to be equal to the the "width" of the loop r w 
that at the same time represents the time that electron 
spend on site zjy+i- Thus, the time dependence of the 
tunneling process is simply given by 



-(£at+i— £o) r n 



(89) 



Making an analytical continuation to real times r w = it w 
and taking the last integral over t w we obtain the delta 
function 



2tt(J(£ 



AT+l 



(90) 



that indeed shows that the process is elastic, i.e. electron 
can tunnel only to the state with exactly the same energy. 

Finally, for the tunneling probability of the elastic pro- 
cess we obtain 



N 



Pel = W £(£a/- + i - 90 Yl P ki 



(91) 



k=i 



where the factor w = n(£o)(l — ^(6v+i) takes into ac- 
count the occupation numbers of the initial and final 
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states. We see that the total probability P e \ contains 
the product of all conductances, mean level spacings and 
inverse energies along the tunneling path. For this rea- 
son it is convenient to introduce the geometrical averages 
of theses quantities along the tunneling path that were 
introduced in Sec. Ill B II Writing the total probability 
P e i in terms of the quantities g, J, E we obtain the result 

CD- 



T N+l 7yv+i 




Ti = -OO TQ ► TQ Tf -OO 



2. Long range Coulomb interaction 

In the presence of the long range Coulomb interac- 
tion the situation is more complicated since the integrals 
over variables r±k cannot be taken on each site k in- 
dependently. However, one can estimate the tunnelling 
probability P e i by finding its upper and lower limits: In- 
deed, neglecting the long range part of the potential (as 
we did above) one gets the upper boundary. The lower 
boundary can be obtained by considering such trajecto- 
ries where all Ar^ are of the same sign. In this case the 
electric field created by the charges on step k does not 
interfere with other charges and thus one can implement 
integrations independently. This gives a lower boundary 
for the tunnelling probability Pk which is about a fac- 
tor 2 smaller than the tunnelling probability Pk obtained 
neglecting the long range part of the Coulomb potential. 
This results in the boundaries of the factor 0.5 < c < 1.0 
in the expression for localization length £ e j in Eq. (|14j). 

The above derivation of tunnelling probability P e i has 
been done for the zero temperatures case. Of course 
the finite temperatures and the presence of phonons are 
necessary for the realization of the Mott-Efros-Shklovskii 
mechanism. Calculation of the tunneling probability at 
zero temperature is justified as long as inelastic tunneling 
processes that are considered in the next section can be 
neglected. 

B. Hopping via inelastic cotunneling 

According to the theory of transport through a single 
quantum dot in the Coulomb blockade regime^, at tem- 
peratures T > y/5Ec the electron transport is dominated 
by inelastic cotunneling processes. In this case the pro- 
cess of charge transfer from one contact to the dot and 
the following process of charge transfer from the dot to 
the other contact is realized by means of different elec- 
trons. As a result, the energy of the out coming electron 
differs from that of the incoming one and the dot en- 
ergy is changed after the each elementary charge transfer 
process. For this reason inelastic cotunneling processes 
are allowed only either at finite temperature or if a finite 
voltage is applied to contacts. 

In order to include the inelastic cotunneling processes 
into the scheme that we described above, first, we have 
to generalize the charge representation to finite temper- 
atures. The main difference is that the imaginary time r 



FIG. 9: Tunneling via inelastic cotunneling between 0— th to 
the N + 1— st grains. 

at finite temperatures is constraint to the interval (0, /?), 
where (3 = 1/T is the inverse temperature. Also the 
expressions for internal energies must be properly modi- 
fied. Using the Matsubara technique one can show that 
the proper expressions for the internal energies can be 
obtained from the corresponding expressions at zero tem- 
perature via the substitution 1/r — > ttT/ sin(7rTr), such 
that, for example, the quadruple internal energy given by 
Eq. (@5J) becomes 

Ul = In [ 2 sin 2 (ttT | At | ) / girT 2 ] . (92 ) 

In general, one should also take into account the static 
charges that formally appear from the integer winding 
numbers in the phase representation at finite tempera- 
ture. However, inclusion of static charges is not needed 
for the description of virtual co-tunneling processes that 
we consider. 

The representation of an inelastic cotunneling process 
in terms of the electron world lines is shown in Fig. 
This process resembles the world line of the elastic co- 
tunneling process shown in Fig. |S| since the electron 
charge is transferred from the site io to the site in+i 
and then back to the site io in essentially the same way. 
However, the important difference is that the diagram 
for the inelastic process is constructed out of quadruple 
loops only. This implies that the charge is transferred by 
means of deferent electrons on each hop. 

As in the case of the elastic cotunneling we assume 
that the sites io and in+i have only one state each with 
energies £o and £/v+i respectively. Similarly to the pre- 
vious subsection consideration, we write down the Green 
functions representing the initial and final states: when 
electron tunnels from site io to site zi, a hole is created at 
site io- Such a process is described by the Green function 
of negative time arguments 

G (t' - T ) = n(&>) e^ T °- T °\ (93) 

The process where electron is located on the site In+i is 
described by the Green function of positive arguments 

G n+1 (t' n+1 -t n+1 ) = [l-nfov+i)] e-^+^+i-^. 

(94) 

The time coordinates and T-k are again confined 
within the temporal interval Ark ^ V^o- This means 
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that the time arguments of the quadruple energies IA!^ in 
Eq. (|92j) are approximately the same and equal to the 
"width" of the loop formed by the world line in Fig. |2l 
t w ~ 7q — r o ~ t at+i — t n+i- Thus, the contributions 



of quadruple internal energies U% to the tunneling prob- 
ability Pi n result in the factor 



sin (ttTt w ) 



N 



(95) 



where g is the geometrical average of the conductances 
over the electron path. Integrations over the intermediate 
times Tfc of the inelastic process can be easily done for the 
case of the diagonal Coulomb interaction Efj = E^Sij. 



1. Short range Coulomb interaction 

In the case of the short range interaction the integra- 
tions over At/c and Ar' k can be performed independently. 
Also the integrations over the Ar^ and Ar^ result in the 
equal contributions given by 



f 

j — ( 



dATk e 



-E£\AT k \+n k A Tk 



2/E k 



(96) 



where the energy Ek is defined in Eq. (JTSJ. Collecting all 
the terms for the probability P^ n one obtains 



Pin (j~w ) 



wg 



2ttT 



2tt jv + 1 \Esm(7rTr w ) 



2N 



-Ar u 



(97) 



Finally, in order to find the tunneling probability P^ n via 
the inelastic cotunneling one has to make the analytical 
continuation in Eq. ([§?)) to the real times r w = it w and 
integrate over t w arriving at 



P — 

J- i n. — 



wg 



AT+l 



+oo 



27^+! 



dt 



2mT 



E sinh(TrTt) 



2N 



-iAt 



(98) 



Here the singularity of the function sinh _2Ar (7rT£) is 
taken in the upper half of the complex plain. At zero 
temperature one can easily calculate the integral over 
the variable t in Eq. ([98)1 by shifting the contour of in- 
tegration in the complex plain either to negative com- 
plex infinity in the case of positive A or to the positive 
complex infinity in the case of negative A. In the first 
case we obtain P^ n = 0. This reflects the fact that the 
real tunneling process with increase of electron's energy 
is forbidden at T = 0. In the latter case, A < 0, the 
zero temperature probability is determined by the pole 
of the function sinh _2Ar (7rTt) that results in the answer 
presented in Ea. ([22)l . 

At finite temperatures the integral in Eq. ([§%)) can be 
expressed in terms of the the Euler Gamma functions 
that leads to the general expression ([!£)). 



2. Long range Coulomb interaction 

In the presence of the long range part of the Coulomb 
interaction the integrations over time variables on each 
step cannot be performed independently. This makes it 
difficult to derive an exact analytic formula for the tun- 
neling probability. However, as in the case of the elastic 
cotunneling, one can easily estimate the upper and lower 
boundaries for this quantity: neglecting the off-diagonal 
part of the Coulomb interaction, as above, we get the up- 
per boundary. The lower boundary can be obtained by 
considering the trajectories where all Ar^ (as well as all 
At£) have the same sign. This consideration constraints 
the values of the coefficient c in the expression for the 
localization length to the interval 1/4 < c < 1. One can 
also see that only the sites that are geometrically close 
to each other can interfere. Indeed, let us consider two 
sites k\ and ki in the tunneling path that are relatively 
far from each other. Interference effects related to these 
sites can be included as a small correction. One can eas- 
ily see that the first order correction vanishes because 
the average electric field on the site k\ created by the 
charges on the site ki is zero and vias versa. Thus the 
contribution of the interference effects falls of with the 
distance as 1/r 2 at least. These insures that there is no 
any dramatic effect related to the long-range part of the 
Coulomb interaction apart form the numerical correction 
in the localization length. The same arguments apply to 
the case of the elastic cotunneling. 



X. MOTT TRANSITION IN A PERIODIC 
GRANULAR ARRAY 

In this section we discuss the metal to insulator tran- 
sition in the regular periodic granular arrays. The Mott 
gap Am for such a system in the nearest neighbor hop- 
ping (quadruple) approximation considered in Sec. I VIII is 
exponentially small as function of conductance, Eq. (|3J), 
but always finite. In this section we consider how this 
result is affected by inclusion of high order tunneling pro- 
cesses that are represented by higher order loops in the 
classical model of Coulomb charges. 

In Sec. IIXI it was shown that at gz <C 1 the contribu- 
tion of higher order tunnellings processes to the tunneling 
provability P at zero temperature is suppressed by the 
factor 



gs 



N 



(99) 



where N is the order of the tunneling process. It is clear 
that the same parameter defines the smallness of the con- 
tribution of the high order loops to the thermodynamic 
potential Q of the classical gas. If the tunneling con- 
ductance is not small, gz > 1, the Coulomb interaction 
of two classical charges is strongly renormalized and in 
this case instead of the bare Coulomb interaction Eq in 
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Eq. ((99)) one has to take the renormalized one - that is 
the Mott gap Am given by Eq. J3J). We come to the 
conclusion that the high order loops can be neglected as 
long as 



A M > g6. 



(100) 



Equation (|100[) can be rewritten in terms of the inter- 
granular tunneling conductance g in the way 



9 <9c 



(101) 



where the critical conductance g c is defined by Eq.Q. 
The samples characterized by the tunneling conductance 
g < g c at zero temperatures are insulators with the Mott 
gap Am given by Eq. The result Q for the critical 
conductance g c coincides with the value of the stability of 
the metallic state obtained from the diagrammatic tech- 
nique from the metallic side of the transition g > g c . 5 
This fact confirms that the value indeed separates 
the metallic and insulating states at T = 0. 

Note that the accuracy of the above considerations is 
not sufficient to establish whether the transition between 
the insulating and metallic states is continuous or dis- 
continuous. The details of the behavior of the Mott gap 
Am in the immediate vicinity of the critical conductance 
g c are not known yet. 

So far we considered the limit of the zero temperature. 
Physics of the Mott transition at finite temperature de- 
pends crucially on whether the transition is continuous 
at T = or it is discontinuous. In the case of continues 
transition at T = 0, at finite temperatures the metallic 
and insulating phases, strictly speaking, cannot be distin- 
guished. Indeed, the Mott transition does not break any 
symmetry and the conductivity at finite temperatures is 
finite even in the "insualting" phase where it acquires 
the activation form, Eq. JSJ). However, in the case of the 
first order phase transition at T = 0, due to arguments 
of physical continuity, the transition has to remain of the 
first order in some interval of finite temperatures. Such 
transition would result in the jump of the physical ob- 
servables even at finite temperatures. 



XI. CONCLUSIONS 

We have investigated Coulomb interaction effects and 
transport properties of granular conductors in the insu- 
lating (dielectric) regime. We have considered both the 
periodic- and irregular arrays and took into account the 
random on-site electrostatic potential. We find the de- 
pendence of the Mott gap on the tunneling conductance 
in the dielectric regime of a regular array. The insulator 
to metal transition takes place at the critical conductance 
given by Eq. (0} where the Mott gap becomes of the order 
of the inverse escape time of the electron from a single 
grain gS. The order of the transition remains, however, 
an open question. For the case of irregular granular ar- 
rays we have derived the Efros-Shklovskii law for the con- 



ductivity of a granular system. Depending on the tem- 
perature the hopping conductivity is dominated by either 
elastic or inelastic cotunneling processes. In the latter 
case the localization length is weakly (logarithmically) 
depends on temperature. In particular, our model and 
the obtained results explain the origin of the ES conduc- 
tivity in the spatially periodic arrays of semiconducting 
quantum dots. 

Finally, we would like to discuss the possibility of the 
observation of the Mott law in granular metals. To re- 
mind, in semiconductors, the Efros-Shklovskii law may 
turn to the Mott behavior with the increase of temper- 
ature. This happens when the typical electron energy e 
involved in a hopping process becomes larger than the 
width of the Coulomb gap A c , i.e. when it falls into the 
flat region of the density of states where Mott behavior 
is expected. To estimate the width of the Coulomb gap 
A c - one compares the ES expression © for the density 
of states with the bare density of states v go , i.e. the DOS 
in the absences of the long-range part of the Coulomb 
interactions obtaining 



A, = 



2d 



(102) 



As we discussed in Sec. [H] in granular metals the den- 
sity of states v go that is relevant for hopping conductiv- 
ity does not coincide with the physical DOS, since only 
the lowest energy excitations in a single grain must be 
counted in determination of v go . 4 

Bare density of ground sates v go can be easily esti- 
mated for the case of strong on-site disorder that we as- 
sumed so far as follows: The width of the potential distri- 
bution jii is approximately given by the average on-site 
Coulomb energy so that 



1/E^a d . 



(103) 



Inserting this value into Eq. (|102|) for the width of the 
Coulomb gap we obtain A c ~ Eq meaning that there 
is no flat region in the density of ground states and, 
thus, the Mott regime is not possible. However, the Mott 
regime may appear in the case of weak electrostatic disor- 
der ( for a detail discussion of this issue see RefA). In this 
case the subgap states will appear very rarely over the 
granular sample and the bare density of ground states 
v go is greatly reduced and as a result the width of the 
Coulomb gap substantially shrinks. 

One can easily repeat all the considerations that led 
us to the expressions for the hopping conductivity for 
the case of flat density of ground states arriving to the 
Mott expression 



a~a exp[-(T M /T)^r], 
where the temperature T^ 1 is given by 



rpM 

1 



(104) 



(105) 
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with the dimensionless localization length £ given by ei- 
ther equation (|14j) or (|21j) depending on the temperature 
regime. 

The crossover temperature T* between the two regimes 
is easily obtained exactly as in the case of semiconduc- 
tors; by comparing Eqs. QJ and (|1()4|) and using the cor- 
responding expressions for T Q we get 



i/(d-i) 



T* - e 2 < 



(106) 



As a final note, we would like to comment on the possi- 
bility of fixing the numerical coefficient in the expression 
for the temperature Tq that enters exponents of the ex- 
pressions for the hopping conductivity. In the case of 
hopping conductivity in semiconductors this coefficient 
can be fixed by means of the mapping the problem of the 
impurity states in semiconductors onto a simple classical 
percolation model that can be investigated completely 
numericall}i2ii22iffi. In granular metals, on the qualita- 
tive level, the situation is very similar, however, it is not 
clear that a model of a granular metal can be reduced 
to a simple unique percolation problem without making 
additional uncontrolled assumptions. The main compli- 
cations are that the tunneling process in granular metals 
(i) in general goes via the curved trajectory, (ii) it de- 
pends on the number of the parameters of the interme- 
diate grains in the tunneling path as E£,5i,gij. For this 
reason, in general, we expect that the coefficient in the 
expressions for To is not universal and is dependent on 
the morphology of a granular sample. 

After the present study was completed, we became 
aware of the preprint where the hopping transport in 
granular metals via elastic and inelastic cotunneling in 
the linear regime ( low applied voltage ) was considered; 
their results essentially agree with ours. 
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APPENDIX A: MINIMIZATION OF THE 
PROBABILITY OF THE INELASTIC 
COTUNNELING 

The result for the hopping conductivity in the regime 
of inelastic cotunneling is obtained via minimization of 
expression ([T8|) with respect to the hopping distance N 
under constraint 



following from the ES expression for the density of states, 
with b being a constant. Minimization procedure can be 
easily implemented presenting the tunneling probability 

Pin 



e , 



(A2) 



with the action A given by 
2A" 



lnr(Ar + iA/27rT) 



+ lnT(N - iA/2irT) - lnT(2N) - A/2T, (A3) 

where & n is the localization length given by (|21|). Taking 
into account that according to Ea. ([Al|) dA/dN = — A/A/" 
minimization of the action A with respect to N in the 
limit of the large N results in the equation 

In 4 + 2 /£ in = ln(l + x 2 ) + nx [1 + (2 /tt) arctan x) , ( A4) 

where the parameter x is defined as 

x = A/2ttTN. (A5) 

We see that Eq. (|A4|) fixes the ratio A/TN for a given 
localization length. Using Ea. (|A4j) and that N ^> 1 the 
action (|A3|) can be simplified to 



A = — 1 H — arctan x ] , 

TV 7T 



(A6) 



that with the help of Eqs. (|A1IA5|) can be further written 
in terms the parameter x as 



A = 



2be 2 
aTR 



xtt [1 + (2 /tt) arctan x) . ( A7) 



Within the logarithmic accuracy we can consider the 
localization length as a small parameter. In this case 
Eq. (|A4|) gives that x « 1 /nt; and from the above equation 
we finally obtain the resulting action 



/ 8be 2 



(A8) 



NaAk/e 2 = b, 



(Al) 



We would like to note that, in principle, one should con- 
sider an alternative minimization scheme for the tunnel- 
ing process: The tunneling first occurs with A = as 
in the case of the elastic process, and then the electron 
gets an extra energy due to phonon or electron collision 
not included directly in the probability of the inelastic 
process under consideration. One c an easily show that 
in this case one would get A = —^/8be 2 /T a^i n k that 
coincides with the action (|ASf> obtained in the leading 
order in the parameter l/£; n . One can check, however, 
that taking large but finite l/& n results in a correction 
to the action (|A8|) that consist of the reduction of the co- 
efficient y/S. Thus the process with intrinsic inelasticity 
that we considered is a bit more favorable. This insures 
that the minimization scheme used by us is the right one. 



21 



1 A. Gerber, A. Milner, G. Deutscher, M. Karpovsky, and 
A. Gladkikh, Phys. Rev. Lett. 78, 4277 (1997). 

2 R. Parthasarathy, Xiao-Min Lin, K. Elteto, T. F. Rosen- 
baum, and H. M. Jaeger, Phys. Rev. Lett. 92, 076801 
(2004). 

3 K. B. Efetov and A. Tschersich, Europhys. Lett. 59, 114, 
(2002); Phys. Rev. B 67, 174205 (2003). 

4 J. Zhang and B. I. Shklovskii, Phys. Rev. B 70, 115317 
(2004). 

5 I. S. Beloborodov, K. B. Efetov, A. V. Lopatin and 
V. M. Vinokur, Phys. Rev. Lett 91, 246801 (2003). 

6 I. S. Beloborodov, A. V. Lopatin, and V. M. Vinokur Phys. 
Rev. B 70, 205120 (2004). 

7 M. V. Feigel'man, A. S. Ioselevich, and M. A. Skvortsov, 
Phys. Rev. Lett 93, 136403 (2004). 

8 V.I.Kozub, V.M.Kozhevin, D.A.Yavsin, S.A.Gurevich, 
Pis'ma v ZhETP, v.80, N5, (2005). 

9 Y. L. Loh, V. Tripathi, M. Turlakov, \c^nd-mat/050174^ 

10 B. L. Altshuler and A. G. Aronov, in Electron-Electron 
Interaction in Disordered Systems, ed. by A. L. Efros and 
M. Pollak, North-Holland, Amsterdam (1985). 

11 L. P. Gorkov, A. L. Larkin and D. E. Khmelnitskii, Pisma 
Zh. Eksp. Teor. Fiz. 30, 248 (1979) [Sov. Phys. JETP Lett. 
30, 228 (1979)]. 

12 P. Sheng, B. Abeles, and Y. Arie, Phys. Rev. Lett. 31, 44 
(1973); B. Abeles, P. Sheng, M. D. Coutts, and Y. Arie, 
Adv. Phys. 24, 407 (1975). 

13 M. Pollak and C. J. Adkins, Phil. Mag. B 65, 855 (1992). 

14 L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975) 

15 I. S. Belobo rodov, A. V. Lopat in, V. M. Vinokur, and 
V. I. Kozub, |cond-mat /0501094 

16 V. Ambegaokar, U. Eckern, and G. Schon, Phys. Rev. Lett. 
48, 1745 (1982); for review see Ref.|27|. 

17 A. Schmid, Phys. Rev. Lett. 51, 1506 (1983). 

18 A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg 
Rev. Mod. Phys. 68, 13-125 (1996). 

19 A.I. Yakimov, A.V. Dvurechenskii, G.M. Min'kov, 
A.A Sherstobitov, A.I. Nikiforov and A. A. Bloshkin, to 
be published in JETP. 

20 T. Tran, Xiao-Min Lin and H. Jaeger, to be published. 

21 B. I. Shklovskii and A. L. Efros, Electronic properties of 
Doped Semiconductors, Spring er-Verlag, New York, 1988). 



22 We consider the system periodic if (i) all the granules are 
of the same size and shape and arranged in a regular struc- 
ture; (ii) all the tunneling conductances between the near- 
est neighbors are equal; and (iii) there is no disorder in the 
dielectric matrix that causes local fluctuations in granules' 
chemical potential. 

23 M. Muller and L.B. Ioffe, Phys. Rev. Lett. 93, 256403 
(2004) 

24 S. Pankov and V. Dobrosavlijevic, Phys. Rev. Lett. 94 
046402 (2005). 

25 Eq. (fTSJ) gives the probaillity of spin unresolved process, 
i.e. it includes the sum over spin components of the final 
state. 

26 N. F. Mott, Adv. Phys. 16, 49 (1967). 

27 G. Schon and A.D. Zaikin, Phys. Rep. 198, 237 (1990). 

28 D. A. Averin and Yu. V. Nazarov, Phys. Rev. Lett. 65, 
2446 (1990). 

29 J. E. Mooij, B. J. van Wees, L. J. Geerligs, M. Peters, R. 
Fazio, and G. Schon, Phys. Rev. Lett. 65, 645648 (1990). 

30 J. S. Meyer, A. Kamenev, and L. I. Glazman, Phys. Rev. 
B 70, 045310 (2004). 

31 A. Altland, L. I. Glazman, and A. Kamenev, Phys. Rev. 
Lett. 92, 026801 (2004). 

32 B. I Shklovskii, Fiz. Tekh. Poluprovodn. (S. -Petersburg) 6, 
2335 (1973) [Sov. Phys. Semicond. 6,1964 (1973)]. 

33 L. D. Landau and I. M. Lifshitz, Quantum mechnics, 3d 
ed., Elsevier Science, Oxford (1977). 

34 We use the dimensionless units for energy of the classical 
model assuming that the effective temperature T = 1. 

35 D. P. Arovas, F. Guinea, C. P. Herrero, and P. San Jos, 
Phys. Rev. B 68, 085306 (2003). 

36 I. S. Beloborodov, K. B. Efetov, A. Altland and 
F. W. J. Hekking, Phys. Rev. B 63, 115109 (2001). 

37 A.A. Abrikosov, L.P. Gor'kov, I.E. Dzyaloshinski, Methods 
of quantum field theory in statistical physics, Dover Publi- 
cations, Inc. New York. 

38 V.Ambegokar, B.I.Halperin, and J.S.Langer, Phys. Rev. B 
4, 2612 (1971). 

39 A. Miller and E. Abrahams, Phys. Rev. 120, 745 (1960 ). 

40 M. V. Feigel'man, A. S. Ioselevich, \cond-mat/0 50248 1\ 



